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O ■ Abstract 



In the context of the p-spin spherical model for generalized spin glasses, we give an estimate of the free 
energy barriers separating an equilibrium state from the metastable states close to it. 



1. Introduction. 

The structure of the phase space in the p-spin spherical model below the dynamical 
transition temperature Td is characterized by the presence of an exponentially high 
number of metastable states [1,2]. If we consider corrections to mean field theory when 
the size iV of the system is finite but large, the dynamics is determined essentially by 
two factors: the mutual disposition of the states and the free energy barriers between 
them. 

The real replica method shed some light into the structure of the metastable states of 
this model [1,3,4] and the same method can be used to estimate the barriers. It is now 
known that, given an equilibrium state, there are many metastable states of various 
energies at finite overlaps with it; it is therefore interesting to evaluate the free energy 
barriers between an equilibrium state and the metastable states close to it. 

In [4] we introduced a three replica potential: the first replica is located into an 
equilibrium state, while the second one is constrained to stay into a metastable state 
at given overlap with the first one (for particular values of the overlap this second state 
can be of equilibrium too); we then calculated the free energy of a third replica as a 
function of its distances from the other two and found two non trivial minima of this 
free energy, corresponding to replica 3 in equilibrium into the state of replica 1 or into 
the state of replica 2. In this context, it is reasonable to give an estimate of the barrier 
between these two states, following the free energy profile of replica 3 when it moves 
from a minimum to the other. This is the purpose of the present letter. 



2. The method. 



The p-spin spherical model is defined by the Hamiltonian 



H(a)= J2 Jii-i^ii ■■■**,, (2- 1 ) 
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where the a are real variables satisfying the spherical constraint jj Yli a 1 = 1 an d the 
couplings Ji lm ..i are Gaussian variables with zero mean and variance 2 jvp-i [5,6,7,8]. 
We consider the case p = 3. 

The organization of the states in this model is quite complex [1,2]. In the range 
of temperature T c < T < Tj (T c is the static transition temperature) the equilibrium 
states coincide with those TAP solutions [9] which optimize the balance between the 
free energy / and the complexity £(/), i.e. solutions which minimize the function 
(/) = f — TE(/); all other TAP solutions correspond to metastable states. 

To inspect the structure of all these states, i.e. their mutual overlaps, we defined in 
[4] a three replica potential Vs(qi2\qi3, Q23), having the following features: replica 1 is 
an equilibrium configuration of the system, while replica 2 is a typical configuration of 
a system forced to equilibrate at overlap qi2 with 1; in other words, replica 2 chooses 
the most convenient configuration compatibly with the imposed constraint; V3 is then 
the free energy of a third replica 3, constrained to have overlaps (qis, #23) with the first 
two quenched replicas. We studied V3 in the plane tc = (913,523), at fixed value of qi2'- 
a minimum of the potential in this plane corresponds to replica 3 having found a state 
into which it can thermalize. 

In the temperature range T c < T < Ta, at given gi2, V3 has two non trivial minima, 
from here on M\ and M2: 

• In Mi replica 3 is located near replica 1, that is in the same equilibrium state; 
indeed, in Mi the free energy and the self overlap of replica 3 satisfy the relation 
of TAP equilibrium solutions. 

• On the other hand, M2 corresponds to replica 3 near replica 2 and its free ener- 
gy and self overlap satisfy the relation of TAP metastable solutions; this means 
that replica 3 has thermalized into a metastable state at distance q\2 from the 
equilibrium state of replica 1. 

Varying 912, M\ always corresponds to the same equilibrium state, while M2 identifies 
with different metastable states, at various distances. 

The important thing is then that, fixed §12, we have the possibility of studying the 
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free energy barrier between two well defined states, via the analysis of the free energy 
contour around Mi and M2 given by V3. In particular, we are interested in the barrier 
that has to be crossed to go from the equilibrium state represented by Mi into the 
metastable state represented by Mi. Indeed, the computation of this barrier is useful 
for a comparison with the dynamical situation at finite N, in which a configuration 
starts at time zero from an equilibrium state and after an exponentially large time 
jumps to a metastable state. 

The plane tt is the image of the phase space T through the mapping which, fixed qi2, 
maps a configuration a of V into its distances (qi3, #23) from configurations 1 and 2. In 
this way, a path in F is mapped into a single path in tt, while, obviously, the opposite 
does not hold. Due to this, a path in the plane n could correspond to no continuous 
path in the phase space and for this reason we can only give a lower bound for the free 
energy barrier between the two states. 

The proposal is to find in n the path linking Mi to Mi which minimizes the variation 
of V3 and to consider this variation as a lower bound for the free energy barrier between 
the two corresponding states. To this aim it is clearly important to consider not only the 
minima of V3 but also its saddles; indeed, as can be seen from next figures, the problem 
of finding the best path from M\ to M2 reduces to single out the chain of saddles which 
minimizes the variation of V3; the estimate of the barrier is then given by 

AF> V 3 (S max )-V 3 (M 1 ) (2.2) 

where S max is the highest saddle crossed along the path. 



3. The results. 

Before analyzing the structure of minima and saddles of V3 at various values of 512, 
we remind that the minimum M2 exists only in the range < q\2 < q(T), and that for 

-4- 



<7i2 — Q*(T), the state associated to M 2 is an equilibrium one [4]. This means that, given 
an equilibrium state, the nearest metastable state that can be seen with this method is 
at overlap q with it and the nearest equilibrium state is at overlap q*. 
We can distinguish three different ranges: 

i- Q12 < Q* : horn Figure 1 we can see that in this range there are only two saddles, 
Si and £2, and the path which links the two minima is then 

Mi -> Si -> S 2 -»■ M 2 . (3.1) 

Since the highest saddle is Si, we have 

AF> V 3 (S 1 )-V 3 (M 1 ) . (3.2) 

It is worth to note that the value of V3 in Mi and Si does not depend on qi2 and 
therefore in this range the barrier is constant. In [3] it has been introduced a two 
replica potential V2, function of the distance qi2 of replica 2 from the equilibrium 
state of replica 1. V2 has two minima corresponding to equilibrium states with 
zero mutual overlap: the most natural hypothesis is then that the maximum sepa- 
rating these minima represents a first estimate of the free energy barrier between 
remote equilibrium states. It turns out that this estimate coincides with (3.2), in 
agreement with the small value of qi2 in this range. 
ii. q* < qi2 < q': at q* S2 becomes higher than Si (see Figure 2) and so now the 
path is the same as in (3.1), but the barrier is 

AF> V 3 (S 2 ;qi2)-V 3 (Mi) (3.3) 

that depends on qi2- We note from Figure 2 that in this range there is a value of 
qi2 at which a new saddle S appears together with a maximum; however, since for 
qi2 < q' we have V 3 {S) > V-j^), the path that minimizes the variation of V3 still 
is (3.1). 
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Figure 1: The contour lines of V3 for <ji2=0.25 and ,3=1.64. The axis are x=^=(qi 3 +q 2 :i) in the 
range [0.25:0.6], and 2/= -7= (913—1723) in the range [—0.3,0.3]. The two minima are visible on the 
right: M\ (up) and M 2 (down); on the left the two saddles: Si (up) and S 2 (down). 
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Figure 2: The contour lines of V 3 for <ji 2 =0.36. The axis are the same as in Figure 1, with 
ranges xG[0.35,0.75] and yE[— 0.2,0.2]. The new saddle S is clearly visible on the right, in the 
middle between the two minima. 

iii. q' < qi2 < q: at q' the path that minimizes the variation of V3 suddenly changes 
because S becomes lower than S% ; so we have 



Mi -»• S -»■ M 2 



(3.4) 



and 



AF> V 3 (S;qi 2 )-V 3 (M 1 ) 



(3.5) 



This situation holds up to q±2 = q, when Mi disappears merging with S. 
The whole behaviour of V3 in the saddles is shown in Figure 3. 
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Figure 3: V3 evaluated in the three saddles Si, S2 and S as a function of <ji 2 . For ,3=1.64, 
g*=0.295 and q'=0.362. 

Finally, Figure 4 shows the behaviour of our estimate of the free energy barrier. 
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Figure 4: The free energy barrier AF as a function of 912- 



4. Conclusions. 

Before analyzing these results it is useful to consider some properties of the potential 
in the stationary points. 

In the computation of V3 we introduced the overlap matrix Q 33 of replica 3, and as- 
sumed for it a one step RSB form, with parameters (x 8 , si, so), where x s is the breaking 
point and s\ > sq are the overlaps. 

The physical meaning of breaking or not the replica symmetry is the usual one adapted 
to this context. An RS form of the overlap matrix can be associated to two very different 
situations: the first one corresponds to a system which finds an exponentially high 
number of states (as in the p-spin model for T c < T < T^), while in the second one 
the phase space consists of just one state (as in the paramagnetic case) . On the other 
side, an RSB form means that the phase space is dominated by a number of order N of 
states (as in the p-spin model for T < T c ). 

It can be shown that in the minima of the potential it holds si = sq, i.e. Q 33 turns 
out to be replica symmetric. This is what we expect: in M2, for example, 523 is big 
enough to constrain replica 3 to see just one state, the one of replica 2, and for this 
reason Q 33 is symmetric, s\ = sq representing the self overlap of the state. 

This symmetry is a very special feature of this kind of situation; indeed, in the 
saddles Si and £2 the matrix Q 33 is actually broken, meaning that the number of states 
accessible to replica 3 is in this case of order N and that replica 3 does not thermalize 
into any of them. 

Surprisingly enough, in the saddle S the matrix Q 33 is symmetric again; moreover, 
in this saddle the free energy fs of replica 3 and the overlap s± = sq satisfy the relation 
between free energy and self overlap of TAP solutions and fulfil the stability condition 
with respect to fluctuations of the overlap (longitudinal stability condition) of [1,8]. 
These facts suggest that S corresponds to a well defined TAP solution, but that, despite 
of the longitudinal stability condition, this solution does not identify a stable state, i.e. 
it is not a minimum of the TAP free energy. We can then make the hypothesis that S 



represents a real saddle in the phase space and that it is just the saddle which links the 
two states corresponding to Mi and M 2 ; this conjecture is supported by the geometrical 
position of S right between Mi and Mi (see Figure 2). If this were true, in the range 
in which the path is Mi — > S — > M2, expression (3.5) would give the exact free energy 
barrier, not only a lower bound. 

We note that the existence among TAP solutions of saddles which satisfy the longi- 
tudinal stability condition of [1,8] requires a more complete analysis of the stability of 
TAP equations for this model. Moreover, it would be interesting to make a computation 
of the complexity of these saddles, to single out their contribution to the dynamics of 
the system. 
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